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ABSTRACT 


The elastic differential scattering cross sections for 


+ + - ~ 
LOD TN 5 8 


aN 5 + Ar were measured for angles of 39° to 


50° with incident lab energy of 100 eV. The apparatus utilized 
the properties of axially symmetric, non-uniform magnetic fields 
(large thin solenoid) to focus all scattered ions with vector 
momentum = within = onto a detector placed down the magnetic 
axis. Thus, the solid angle was increased by a factor of about 
300 over conventional methods with comparable resolution in AQ 
(about 1°). A typical value for @ (43°) lomo,2 x ire” ene 

H 


+ 
The H data was compared with data of Magnuson, et.al., viaa 


scaling law and was in excellent agreement. A typical value’ 
oa -20 2 : : 

for 0 (43 ) Ss 25a cm. A Simple refinement of a 

H 

2 

+ 
classical model for H, scattering by Feist relates these cross 
ti by @ (0) = PcG (0) = Be DAG fone here P is an 

sections by Fe ) ( = r where 


H5 Hie nL B. H 


orientation factor, D the binding energy, and £ the internuclear 
separation. This expression relates the cross sections to within 
the accuracy of the experiment thus verifying the classical 


approximations. 
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I.  TNTRODUCT ION 


Atomic processes such as the scattering of complex ions 
by heavy atoms are properly described by quantum mechanics. 
Unfortunately, for a many-body problem the mathematics becomes 
untractable. In order to produce useful solutions approximations 
must be introduced. One class of approximations involves the 
use of semi-classical models of the collision process. By 
Sacrificing rigor, useful solutions are obtained that still 
provide insight into the physical process involved. 


The classical model investigated here is used to calculate 


+ 


the differential scattering cross section of the Hy 


10on by a 
spherically symmetric potential. The theory was developed by 
Feist [1] based on the earlier work of Bates, Cook, and 

Smith [2] as follows. First, calculate tHe differential 
scattering cross section of a single proton from a spherically 
Symmetric potential. Then calculate the probability that the 


second proton is oriented with respect to the first such that 


a a e ° 
the H, molecule is not disrupted by the collision. The dif- 


+b 


> 1S the product of 


ferential scattering cross section for H 
this orientation probability times the cross section for the 
single proton. 

The purpose of the experiment reported here was to test 


the theory as developed by Feist. This was done by measuring 


: + 
the differential scattering cross section of H by Argon and 


a0 by Argon at incident energies of 100 eV and laboratory angles 


2 
from ayo) tae 49°, Thus, the orientation probability was deter- 


mined experimentally and compared with Feist's calculation. 


10 





Il. THBORY 


A. THE MODEL 

This classical model is applicable to the scattering of the 
a ‘ ; : , 
H,, 10n Over a range of incident energies Ey (in laboratory 
coordinates) by a spherically symmetric potential such as heavy 

+ + 

noble gas atoms. The reaction Studied was H, ; At = A, + Ar 
where the underline indicates the fast particle in laboratory 


e 


coordinates. The two protons of the H, 


ion have an equilibrium 
Separation in the ground state of 1.06 angstroms and a disasso- 
ciation energy of 2.648 eV. The lower energy limit where this 
model could no longer be expected to apply would be where the 
incident ion energy becomes comparable with the ion binding 
energy. This requires that E,>> 2.648 eV or B. > 30 eV. The 
upper energy limit where this model could no longer be expected 


to apply is governed by the uncertainty principle. In this 


application it is required that the de Broglie wavelength be 


much less than the distance of closest approach or rt o7\, = — 
oO 
which may be written as r oe - Now r_ may be further 
T2mE, i 


related to Eo by considering a head-on collision with the shielded 
coulomb potential for Argon (this potential is discussed in Section 
II B, and Appendix A). Considering only the dominant term of the 


potential one can write 


Where A and @ are parameters of the potential. Substituting 


1 


values for the conStantsS one obtains 


l 
2 
1g 
O 
" = Lt eee 
and Eel) — oe - 
O ac 


with E, in units @f eVvgand t in units: of Bohr radiigee fhe value 


l, 


OF EY that satisfies these equations is Breen 1. 13tehO.-cVeeeDhms 


shows that a classical model for the scattering of H, should be 


valid for incident energies up to about 10 KeV. 


+ 


In a scattering event the H, 


10n may be disrupted due to the 
following two effects. First, consider the two protons approach- 
ing a spherically symmetric potential with the same impact para- 
meter, b. Using the concept of an impulse approximation, con- 
sider each proton-target collision separately, yielding a de- 
scription in terms of two two-body collisions. Due to the 
spherical symmetry each proton scattered is confined to a plane 
defined by the scattering axis and a line from the scattering 
center to the proton. The angular difference between the 
scattering planes, A® will result in a difference in velocities 
of the two protons after scattering such that their relative 


‘ : a Nas 
energy may be greater than the binding energy and the H, ion 


Zz 
will be disrupted. Figure 1 illustrates this process. 


Second, consider the protons oriented such that they will be 


scattered in the same plane but with different impact parameters. 


Ik2 


a ee 
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This will result in a difference in scattering angle AQ and 
therefore again a relative velocity between the two protons. 
The relative energy may again exceed the binding energy so that 
the ion will be disrupted. Figure 2 illustrates this process. 
In the actual case these two processes are not separable 
and disruption may occur due to the combined effect of AY and 
A®. The necessary condition for the ion remaining intact after 
collision is D> p(Aav)* where D is the disassociation energy, 
b= “e the reduced mass of the two protons, and AV is the 
difference in velocities due to both AMY and A®. Whether the 


H, war = H, 2 £m Yeaction will o€cur, is thenya tumgotion ehethe 
orientation of the two protons before (or equivalently after) 
the collision. The differential scattering cross section for 
the H, ne). oe may be calculated as follows. For E>>D 
treat the cet fete of the H, ion as two events consisting of 
the scattering of proton "A" and proton "B" (impulse appr oxi- 
mation) ._ Calculate the differential cross section for proton 
"A", © (8). Calculate the frobabilagey,, I, othaitkeproton es e1- 

‘ 2 
ordented = with méspect to proton "A" suchéBhat D> 3 U( AV) st 
then follows that | 


Q) = OM. A 
atl Vier o £9) (1) 


The electron is not treated as a scattered particle in this 
model but may be thought of as simply providing the binding 


energy. 


14 
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B. CALCULATION OF 90 ,(9) 
H 


The galculatien er oa was done on an IBM 360-67 digital 
computer using a program written by Feist P22 tne program 
utilizes standard numerical techniques to solve the classical 
differential scattering equations. The potentials used are 
shielded coulomb potentials from Byatt [31 of the form 

Z -b x 
(2 ce ) 


n 


ze 
Yr 





V(r) = 


with , Cc, = 1 but Cand db. otherwise adjustable. Byatt lists 

two potentials for Ar a@rmone with and one without correction 

for polarization. The two potentials yielded cross sections” 

for H” that were within 13% of each other for an incident energy 

of 100 eV and scattering angles from BO to som (laboratory co- 

ordinates). Due to the close agreement of the two potentials, 

only the one corrected for polarization is presented here. 

Byatt indicates the potentials should be accurate to within 5% 
-b r 

for 3 Gace  >.0i%. ‘For the region where the potentials were 


n 
_ =b © 


used it Mas notedethat & SF e 7: >.052. Thus it was felt that 
n 


these potentials should be satisfactory for calculating 0 (9). 
H 


C. CALCULATION OF THE ORIENTATION PROBABILITY 


oar oe ; 
The energy criterion for Bihe H, 10n remaining intact may be 
put as M 
a 2. Pp é 
. 2 CAME abe? 4 ee) 
1. 
ox AV 2 (2) 
max M 
Pp 
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Considering the case where the target mass 1S much greater than 
oe cee _ 
the mass of the H, ion, the H,, i10n has energy after collision 


2 
ES = me Nes both in the laboratory and center of mass systems. 


Substituting this into equation (2) gives 


1 
AV a Ce ie 


max E oO 
Oo 


Now it was shown in Section II A, that OW cx ® file, Ab, Ag). 
Therefore, for a fixed incident energy and impact parameter, 
there exists a function g such that A® = g(Ab). This defines 

a closed curve ina plane perpendicular to the scattering axis 
that proton "B" must pass through in order for the ion to remain 
intact. The area enclosed by the A® = g(Ab) curve is designated 
py S. 

The axis of the H, ion 1S randomly oriented with proton ‘'B" 
having no preferred position with respect to proton "A", The 
location of proton "B" will therefore define a sphere about pro- 
ton "A'' of radius £, the internuclear separation. The orientation 
probability is then given by the ratio of § to one-half the area 


of this sphere 





= ; (&) 


The justification for using the area of a hemisphere in the 
denominator of equation (3) is that the "other half" is re- 
dundant. No exact solution has been found for equation (3) and 


the following approximations were made. 


1 


First, it has been assumed that the target is infinitely 
— : 
heavy so that the speed of the H, 10n 1S unchanged by~the col- 
lision. For the reaction studied experimentally, this is valid 


in that M , = 2.032 amu and M,. = 39.948 amu. For example, a 
H 


fh 


Zz 
100 eV a 10n scattered at dan? has an initial speed of 9.80 x 10 


i 


m/sec and a final speed of 9.75 x 10° m/sec. It is further 
assumed that the two protons will both have this same speed after 
collision and that disruption will occur due to difference in 
direction only. Feist Li] showed that the disruptive effect 

of the difference in the magnitudes of velocities of the two 
protons is much less than the disruptive effect of their dif- 
ference in direction. For the case where A@ = 0, and Mest 


A@ 1s small, AV =v [lao | Using equation (2) this 
max O 


max” AS = O° 
becomes 
_ f4n]% 
LAO dag =O Ei ’ (4) 
Figure 3a illustrates this case. For the case where AQ = O 


(or equivalently Ab = O) and assuming A@ is small, 


ae = Va Sin 9 Uae ds — Using equation (2) this 
becomes 
qj _ | 4D ]% 1 
Ce tela = Te ze cane (5) 


Figure 3b illustrates this case. Using the classical definition 


fox differential scattering |2nbdbl = lo (9) 2 on Crd", 10 
H 
1s seen that o (9) sin QO 
yt 
i a Q 
AD ax! ag = 0 b [4 ae} AG = 0 ” 
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scattering plane. 


Lo 


which, upon combining with equation (4) becomes 


,(9) sin 0 
H ds 
| Ab ax! ag 0 a ay + — 


Equations (5) and (6) give maximum values for either Ab or AQ 


a : : 
such that the H, ion is not disrupted. The 4% = g(Ab) curve 


must therefore be enclosed by these maximum values. An upper 
bound on the area, S,* enclosed by the AM = g(Ab) curve is there- 


fore 


i]. 


S =p 
maxX m 


JL2b( ag 


athe = 0 ewan =o 


Substituting equations (5) and (6) into the above and using 
equation (3) an upper bound on the orientation probability is 


found to be 


Co a) 
P = § D em | 

max BG a (7) 

and 0 (9) = ig HO Ie (8) 
H»5 Tikes mo B He 


Due to the simple form of equation (8), it is calculated numeri- 

cally in the same computer program used to calculate O (9). jet 
H 

should be further noted that D is the disassociation energy for 

ground state and £ the corresponding equilibrium distance. Ex- 

cited electronic states for the H, molecule are eee unboowe 

or have very shallow minima in their potential curves. Therefore, 


if the H, 10n 1S excited to higher electronic states during the 


collision process it is almost certain to be disrupted due to 


20 


I ,llsi—(isti—Cser—SNCisiCri 


: , — eames : 
velocity differences after the collision. If the H, 10n is 
“excited to higher vibrational levels, then D would be decreased 
and £ increased. The molecule could still remain intact but 
would have a smaller value of P and hence 0 ,(9)- 

ni 
D. COMPUTER CALCULATION OF 0 oer) AND 0 are’) 
H H 
2 
Figure 4 displays the results of the computer calculations 
for the differential scattering cross sections. One set of 
+ + , 
curves is for the H_ + Ar ~ H_ + Ar reaction at incident energies 
of 100,300, and 500 electron-volts. The second set of curves 
+ 


a : 3 : 
is for the He twat. H, + Ar reaction also at incident energies 


of 100, 300, and 500 electron-volts. The H, curves were com- 
puted using equation (8) and represent a predicted upper limit 
on the scattering cross section. 

It is noted that all of the curves are monotonically de- 


creasing functions of scattering angle and reflect no atomic 


strueture. 


Zi 
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Figure 4. Log of scattering cross section vs. scattering angle for 


H*+ Ar—> H* +Arand H,” + Ar—+H, + Ar. 
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Ill. EXPERIMENTAL APPARATUS 


A. GENERAL DESCRIPTION, 


The experimental apparatus used to measure the cross sections 


+ + ~ + 
fom the H , H, + AX ~ Ou 


3 + Ar reactions is discussed in detail 


by Bush L4]. An overall schematic of the apparatus is shown in 


Fig. 5. The operation of the apparatus was as follows. Three 


- 


“- 
species of hydrogen ions (H , H, > 


Hy) were produced in the duo- 
plasmatron i10n source and were accelerated electrostatically 

to the desired energy. The ions entered a massS analyzer where 
the desired 10n species was bent through a hs? angle to travel 
down the geometric axis of the apparatus. The selected beam was 
focused and entered the scattering cell where the target gas, 

in this case Argon, was contained. The scattering cell was 
centered in a large focusing magnet that produced a known non- 
uniform but axially symmetric magnetic field. The trajectory of. 
an ion with a given vector momentum could then be calculated and 
the detector was appropriately placed to collect this ion. Due 
to the finite detector size, ions scattered within a finite AQ 
(on the order of i) were also detected. It is further noted 
that essentially all the ions scattered ein angle 9 were col- 
lected regardless of their azimuthal angle of departure. Cross 
section measurements at different scattering angles could be made 
in two ways. Either the detector could be fixed on the magnetic/ 


geometric axis and the focusing magnet current varied to select 
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the desired scattering angles, or the focusing magnet current 
could be fixed and the detector moved along the magnetic/geometric 
axis to select the desired scattering angles. In this experiment, 
the latter method was used. The detector system consisted of a 
magnetic electron multiplier whose output was fed into an elec- 
trometer. The entire system was maintained at a vacuum on the 


5 


order of 10 ~ torr during operation by two diffusion pumps. 


B. THE DUOPLASMATRON ION SOURCE. 

The duoplasmatron 10n source was constructed and described 
in detail by Carter [5]. Figures 6 and 7 display the construction 
of the source and the circuit diagram for its operation. The 
Operation of the source basically consisted of maintaining the 
plasma discharge in the Z-electrode cavity and extracting ions 
at the desired energy by placing this potential on the anode. 
Focusing the extracted beam Was accomplished by an electro-magnet 
surrounding the Z-electrode cavity and a conventional Einzel lens 
and deflector stack. 

To increase the beam intensity an acceleration-deceleration 
system was utilized. This consisted of extending the extraction 


lens, L to enclose the entire focusing stack. This enclosure 


deg 
was called the "sock". (See Fig. 7.) By floating the sock at a 
negative potential V (typically-500 volts) and by increasing the 


potentials of the interior lenses and deflectors by the same 


amount, the extracted ions transversed this region with an energy 
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(E + V)e, yielding a higher space charge limit for the beam 
density. Upon exiting from the sock, the hydrogen ions were 
decelerated to energy Ee. This arrangement increased the beam 
intensity at the Scatt@inpgece!] by a factor Gheabout 40. 

A typical maximum beam intensity for mass analyzed protons 
of 100 eV delivered to the scattering cell with 0.5 cm entrance 
aperture at 80 cm from the source was 1 x 1077 amperes. For 
H, ions under the same conditions the maximum beam intensity was 


on the order of 2 x 1077 amperes. The energy spread of the beam 


was estimated to be aboutweemwev for a 100 eV beams: 


C. THE MASS SPECTROMETER. 

The mass spectrometer used to select the desired ion was 
constructed, calibrated and described in detail by Strohsahl L6]. 
Figure 8 is a schematic of its arrangement in the overall appa- 
ratus. The mass spectrometer utilized an electromagnet to pro- 
duce a magnetic field of up to 4 kilogauss confined between two 
4o° Sag e- aoe pole pieces. This was sufficient to mass 
analyze atomic masses up to 40 with energies up to 1 KeV. The 
resolution aa was about 40. A bending angle of 45° was 
utilized and this altered the selected ion beam trajectory from 
the axis of the duoplasmatron to the geometric axis of the main 
scattering apparatus. Operation of the mass analyzer consisted 
of merely of setting the correct magnet current to focus the 


desired ion into the scattering cell. The pole faces and region 
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where undesired species were "dumped" were coated with aquadag 
to prevent static charges from building up. Even so, this region 
required thorough cleaning after about 20 hours of duoplasmatron- 


analyzer operations. 


D. | THERESCATTERING CELL. 

The scattering cell was designed by Bush L4]. Figure 9 shows 
a schematic of its construction. The functions of the scattering 
cell are to (1) allow the collimated incident ion beam to pass 
through the scattering region and be measured by the beam collector, 
(2) contain the target gas in the scattering region and measure 
its pressure, and (3) allow scattered ions to escape over as large 
a range of scattering angle and azimuthal angle as possible. The 
design of the scattering cell was a compromise between these re- 
quirements. 

The incident beam enters the scattering cell through an 
adjustable aperture (diameter from 0.1 to 1.0 cm) which was 
usually set at 0.5 cm. The beam passed through the front section 
which contains the scattering region. Now the great majority of 
the incident ions will not be scattered at all and will continue 
on into the rear section where they are collected either by the 
rear section wall or the beam collector. The incident beam was 
measured by two Keithly Model 410 Micro-Microammeters. The beam 
collector, which collects about 90% of the beam, had three grids 
that prevented slow ions from being collected and suppressed 
secondary electron emission. The beam collector assembly could 


be swung off axis allowing a clear path through the scattering cell 


for system alignment and detector calibration. 
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The target gas was bled into the scattering region through 
a Granville-Phillips variable leak valve. The pressure was 
measured by a VG-1A ion gauge manufactured by Vacuum Industries. 
Since operation of the i0n gauge interfered with the incident 
beam measurement, the pressure was measured before and after 
each run and averaged. Vacuum Industries lists the correction 
factor“fior Argon (relatiVewto dry air) as 0.617. That is, 
Pao 3 0-01 / = P indicated’ 

The rear section of the scattering cell was adjustable along 
the axis of the apparatus and thus regulated the target thickness 
between 0.1 and 0.8 cm. This was set at 0.3 cm for the present 
experiment. 

The scattering angle limits were defined by the slope of the 


, Oo 
front and rear section shoulders of 49° and 36° respectively. 


This range of scattering angles could be extended somewhat by 


considering an effective target thickness. Bush L4] describes 
this process. Ideally, ions scattered into all 360° of azimuthal 
angle could be collected at the detector. However, interference 


by mechanical support required for the scattering cell reduced 


the available azimuthal angle to 290°. 


E. THE FOCUSING MAGNET. 
The focusing magnet consisted of five concentric spools of 


copper wire in a parallel-series arrangement. Each spool was 
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water-cooled by a \ inch copper tubing jacket. A Sorenson 
Nobatron DCR 150-15 power supply maintained the desired current 
setting for periods of several hours with essentially no drift. 

The magnetic field on the center axis was carefully surveyed 
by Kelly L7] using a Hall probe. Figure 10 shows the result of 
these measurements. Kelly also showed that the magnetic field 
was Symmetric about the coil axis to within 0.1%. This allowed 
the axis of the apparatus to be aligned with the coil axis so 
that the magnetic axis coincided with the geometric apparatus 
axis. 

In order to better understand the function of the focusing 
magnet, it is instructive to consider the simpler case of a 
uniform axially symmetric magnetic field. This would correspond 
to the entire apparatus being enclosed in a long solenoid. For 
an 10n of charge q with momentum p and scattered at an angle 9 
the trajectory is a helix such that after one orbital revolution 
the 10n will again intercept the axis at some distance Zo° 
Using these boundary conditions the Lorenz force equation 
F = q v x B can be integrated to give Zo = a Dp €s 9. A de- 
ieetor placed at Z, will then detect the scattered ion. (See 
Fig. 11.) This equation also shows the two ways of detecting 
Scattered ions over a range of 90 and p. Either maintain B con- 
stant and vary the detector position or keep the detector fixed 


and vary B. 
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Figure 10. Magnetic Field Strength at Various Distances 
Along the Magnetic Axis for Various Currents 
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Now in order to measure a differential scattering cross 
section, the solid angle that the finite size detector "sees" 
must be known. For the axially symmetric case this is 
dQ = 27 sin 09 dO. To evaluate d@ consider the various para- 
meters affecting Z 6: the axis intercept distance. A functional 


relationship exists such that: 


Za = £(E,9,B,P,7,C) 
where E = enerom® of incident 10n 
QO = angle of scatter 


B = magnetic field 


radial target displacement from the system 
ax1s 


Re) 
I 


displacement of target atom from the 
Scattering center along the system axis 


3 
II 


C = constant for each reaction including mass 
ratios, inelastic energy losses, etc. 


Differentiating this expression gives 

A 7 —_ 5 AE + Be a TO| is 4| ae A [P20 A Wale 
o *\ OE 50 3B api} -yaltey jim lac 

For a specific measurement AB and AC are zero. Since analytic 

expressions for the various derivatives are not known, they must 

be evaluated numerically by varying E, 9, P, and 7}. Bush, Cook, 

Heinz, and Rodeback [8] estimate that 4p and 47 which are due 

to the finite volume of the scattering region affect AQ by no 

more than 10%. Approximating Ap and An as zero the above ex- 


; | a0 OZ, 
pression can be rearranged to AQ = 52, | AZ, - SE AE 
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The effect of a small (2%) energy spread was studied and it was 
concluded that the "effective" AQ was not changed. The angular 
resolution of the apparatus was degraded somewhat in that not 

all particles scattered within AQ of 9 were collected but this 
was compensated by the fact that some particles scattered outside 
of AQ were collected. The solid angle can therefore be approxi- 
mated as 


woe aeuae) 
aa ~ an sino [(35-) 1 az, 


AZ 





where + is the maximum deviation of the trajectory inter - 
cept from Z, that will still allow the particle to be collected 

by the detector. AZ. can be calculated by considering the de- 
tector geometry and the angle C.. with which the particle inter- 
cepts the Z-axis. As an example of how the approximate expression 
for the solid angle is applied if an analytic expression for 

Zo = £(9) is known, consider again the uniform axially symmetric 


magnetic field. Recalling that Z = = cos 90, it follows that 
aQ~ Baz. 

— p = 

For the non-uniform axially symmetric field, the procedure 
for computing trajectories and solid angles was essentially the 
same as for the uniform case except that the calculations were 
done numerically on an IBM 360-67 digital computer. The program 
that accomplished this is explained in Appendix A. A very brief 


outline of the program is as follows. 


5% 


(1) Express the magnetic field along the axis as a 12th 
order polinomial. This gives Be > (1) Boz” where 
I 1s the magnet pe ee 

(ZamUsSing @this ool aries a vector potential is found. 

(3) Combining the vector potential with the Lorentz force 
equation gives a ond order non-linear differential 
equation. This is expressed as two simultaneous, first- 
order differential equations which are integrated by 
the DHPGC Subroutine of the Scientific Subroutine 
Package for the IBM 360-67 Computer. the trajectory, 
Zo? and OF are now known. 

(4) The above procedure is done over a selected range of 
sattering angles anda olen expressing 9 = £(Z.) 
is generated. 

(5) From 9. and the detector geometry, AZ 1s calculated. 

(6) The solid angle can now be calculated as 
dQ (6) = an 5) sin ate \Z,- The factor = 
1s present due to the construction of the scattering 

Cell... 
Figure 12 shows a series of trajectories computed for H” 


ions with incident energy of 100 eV elastically scattered by 


Argon. 
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FO eta | DETECTOR, 

The detector system employed in this experiment consisted of 
a Bendix Model 306 Magnetic Electron Multiplier (MEM) which was 
fed into a Keithly Model 640 Electrometer. Figures 13, 14 and 
15 show the detector geometry, a schematic of operation, and 
a circuit diagram. ThewgrineipBevol operationsefeune multiplies 
was as follows. An incident ion was accelerated to about 1500 ev 
and ejected an electron from the tungsten cathode upon impact. 

A crossed electric and magnetic field (alnico permanent magnets) 
imparted a cyclical trajectory to the electron(s) down the dynode 
strip. The glass dynode and field strips had a high resistance 
coating @f tin oxide and antimomyl that sustadned the electric 
field. Hach time the electrons struck the dynode strip multi- 
plications occurred so that when the electrons were collected by 
the anode, a gain of about 10+ had been realized. 

Unfortunately the MEM had to be operated in the focusing 
magnetic field. This resulted in a sideways drift of the elec- 
trons moving down the dynode strip so that all of the electrons 
generated were not collected at the anode. This in fact made 
the MEM gain a very Sensitive function of its position in the 
focusing magnetic field. The procedure used was to calibrate 
the MEM under the same environment that existed —e the actual 


cross section measurements. Specifically, this consisted of: 


+ 


3 ions with a focusing magnet cur- 


(1) measuring the gain using H 


Wy 


rent of 6.6 amperes with the detector position, Zo? varied from 
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Figure 13. Dector Aperature Geometry. 
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Figure 14. The bendix model 306 magnetic electron multiplier 
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Figure |5. Circuit Diagram for Electron Multiplier. 
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30.5 to 47.5 centimeters from the scattering cell, (2) measuring 
the gain using H’ ions with a focusing Magmetwcurrent of 46 
amperes with the detector position, Zo? Varied from 26.5 to 427 
centimeters from the scattering cell. Since the focusing magnetic 


field was a known function of current and position, the two gains 


oF 


are plotted against magnetic field in Fig. 16. The gain for H, 


10nsS was measured before and after the cross section measurements 
and was found to remain constant. The gain for H” 10ns was 
measured only once. 

Due to the steep slopes of the gain curves, a small error 
in setting detector position or focusing magnet current will 
result in a large error in the detector current reading. This, 
coupled with the tedious task of calibrating the detector for 
every Species, energy, and magnetic current that a cross section 
ber cones! 1s Fo be nade made this detector system less than 
ideal. For further experimental work with this apparatus, a 
different detector system is being considered. This may consist 
of counting the output pulses from the MEM rather than integrating 


them with an electrometer. 


Ll 
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Figure I6. Multiplier gain vs. magnetic field. 
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IV. EXPERIMENTAL RESULTS 


A. LABORATORY CROSS SECTION MEASUREMENT. 
The equation used for calculating a differential cross section 


in the laboratory was 
I_/G 
Wa 


do(Q) = cm 
( ) my Tasaie x P x Comm hae t or) x t 0d (9) ss 
L Sc 


where I. = detector current measured by an electrometer 


G = MEM gain 
I. = Incident ion beam. Sum o@ the currents to 
Beam Collector and Rear of Scattering cell. 
16 5 ; 
p= 3.54 x 10 particles/cm~-torr based on ideal 
gas law. 
P = Pressure in scattering cell as measured by 


sc 


VG-1A i10n gauge 


Corr. Factor .017; corrects ion gauge reading for Argon 
t = target thickness. Normally equal to .3 cm 
in this experiment 
d91(Q9) = Solid angle computed by digital computer. 


a : A P 
A sample calculation for the H + Ar reaction showing typical 


experimental values is 


~9 P 
(CHR) Pace sa We SAT oes 10 . 


H (lmecpie x 107° amp.)(3.54 x VO neenme dames | 


1 
fleS x 10°torr x -617)(0.3 cm)(0.109 sterad) 


eo? x Moo eae 
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B. MEASUREMENT OF THE H + Ar ~ H + Ar CROSS SECTION. 

THe differential scattering cross Section Mor this remetion 
was measured from 202 to 50°. Five runs were made and the 
data points are displayed in Fig. 17 along with the computed 
cross section using the Byatt potential. Figure 18 shows the 
average of the 5 runs along with the computed cross section. 
Important parameters for the measurements included: 

(1) Incident energy of H ions equal to 100 eV. 
(2) Focusing magnet current equal to 4.8 amps. 
(3) Solid angle typically equal to O.1 steradians. 


(4) Argon pressure in scattering cell typically equal to 
1 x 10°? torr. 


Figure 18 shows that the experimental curve, while having 
about the same overall slope as the computed curves, differs in 
magnitude by*a factor of about 6. "This difference, while not 
catastrophic, was greater than expected. First, there may have 
been a systematic error in the experimental apparatus. A pos- 
sible systematic error could have been an error in measuring the 
MEM gain. The gain was measured with an incident beam several 
orders of magnitude greater than the beam of scattered particles. 
Thus, if the MEM gain were a (slowly) increasing function of input 
current, the measured gain would be greater than the gain for the 
smaller beam of scattered particles. This would result in the 


measured cross sections being too small. 
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Figure 17. Experimental measurements of H*+Ar—+ H* + Ar at 
lOO eV. 
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Figure 18. Average of experimental measurements of 
H* + Ar—H*+ Ar at 100eV. 
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Second, the potential may have been inadequate for this 
application. Varying the potential parameters revealed that the 
computed cross sections were relatively insensitive to small 
changes in the potential parameters. For example, varying the 
dominate terms in the potential by 30% such that the collision 
was made softer decreased the computed cross sections by a factor 
of only about 2. It was felt that the potential parameters were 
at least this accurate. It is noted, however, that the potential 
is Hartree-like and therefore will reflect no atomic structure 
or inelastic process. At the large angles where this experiment 
was conducted, the collisions are relatively hard so that the 
Argon atom may well be excited. The ny 10on would then have .- 
Significantly less energy and momentum after the collision com- 
pared to the elastic case. Since the focusing magnetic field 
acts as a-momentum analyzer, only elastically scattered H” 1ons 
were collected by the detector. Thus, the presence of an 
inelastic process could decrease the elastic cross section by a 
factor of 6 and not be accounted for in the Hartree-like potential. 

The work of Aberth St Lorents [11] tends to Support the 
second analysis. Aberth and Lorents measured the scattered signal 
profile of He Oneonta 10° and 400 eV with an energy analyzer. 
The result shows that for He On Ar, an inelastic event is more 
likely "than elast?eggecattering+by a factompof about 2%. The Gmergy 


um : : 
lost by He in inelastic events was on the order of 10 to 15 eV. 
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Aeplying this to the caselot H- on Ar, 1t seems reasonable to 
expect protons to be a more efficient projectile for exciting 
Ar than He and therefore could account for the factor of 6. 
Furthermore, inelastically scatter protons with energy after 
collision of 85 - 90 eV would not be collected by the detector. 
The choice between the computed and experimental cross 
sections was accomplished by comparing the data in Fig. 18 with 
data supplied by Magnuson, Carlston, and Snyder Pov Magnuson 
ett also measuredstne cross Seclichpof on Atyalthengh stheir 
measurements were made for scattering angles less than 20° and 
energies ranging from 150 eV to 800 eV. Meaningful comparison 
of the data was possible by using a reduced plot outlined by 
Smith, Marchi, Aberth, Lorents, and Heinz L10]. This reduced 
plot consists of plotting T = EO versus 9 = Osin 9 0(9). By 
the use of this scaling law, cross section curves obtained at 
various energies are reduced to a single line. The reason is 
that the product EO represents essentially the same encounter 
with the atomic potential. By comparing the data in this manner 
it was found that the two experimental curves agreed to within 
a few percent and the unaveraged data infact overlapped. Thus, 
it was concluded that the experimentally determined cross sections 
were probably correct. It is further noted that the data of 
Magnuson, et.al. covered a much larger range of T = EO than this 
experiment and does reflect atomic structure such as discussed 


in the preceeding two paragraphs. 


Silt 


C. MEASUREMENT OF THE H, + Ar 7 Hp + AY CROSS SECTION: 
The differential scattering cross section for this reaction 

was measured from econo to 49.7°. Five runs were made and®the 

data points are displayed in Fig. 19. The average of the 5 runs 


1s displayed in Fig. 20. Important parameters for this measure- 


ment included: 


+ 


(1) Incident energy of H, 


10n beam equal to 100 ey. 

(2) Focusing magnet current equal to 6.6 amperes. 

(3) Solid Angle typically equal to 0.1 steradians. 

(4) Argon pressure in the scattering cell typically equal to 


| 5% 1072 


eOLik. 

From Fig. 20 it is seen that the overall slope of the experi- 
mental curve is essentially flat. That the cross section does 
not decrease significantly with increasing angle is possibly due 
to the following instrumental effect. Some particles scattered 
outside of the computed solid angle may reach the detector through 
multiple collisions with the walls of the apparatus and still have 
sufficient energy to overcome the detector grid bias. The sige 
nal due to these extraneous particles could not be subtracted 


out in that they essentially behaved like the "real" particles 


scattered into d8. That is, the signal due to extraneous 


+ 


2 beam intensity and 


particles was proportional to the incident H 
to the pressure in the scattering cell. Normally, this extraneous 


Signal is negligable compared to the signal from particles scat- 


tered within dl. However, in this case at the larger angles the 
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Figure I9. Experimental measurements of He + Ar-H Ar gt lOO eV. 
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Figure 20. Average of experimental measurements of He +Ar+H? + Ar 
at 100 eV. 
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detector signal was quite small, typically on the order of 

-12 
Oo X oO amperes. Therefore, perhaps up to 50 percent of 
the measured signal could have been due to these extraneous 
particles. This would also account for the wider dispersion of 


data at the larger angles. 


D. EXPERIMENTAL RESULTS AND THE CLASSICAL MODEL. 

Since it could not be assertained with certaintly as to 
whether the computed or experimental cross section for Hq" on 
Argon was more accurate, the analysis of the classical model 
was done from both points of view. 

First, it was assumed that the computed cross section for H” 
was correct. The difference between the experimental and com- 
puted cross section was then attributed to a systematic error 
and the experimental H data was normalized to fit the computed 
curves. This was done by moving the H” data wp by A factor of ©. 
Now the H, on Argon cross section measurement had the following 
Similarities with the H” measurement: 

(1) Approximately the same focusing magnetic field (see Fig. 16). 
(2) Approximately the same excursion in detector position. 
(3) Approximately the same calculated solid angle. 
(4) The gain curves for both measurements have about the same 
shape (see Fig. 16). 
It was therefore concluded that if a systematic error were pre- 


+ 
sent, it should have affected the H, measurement by about the 
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‘ + : 
same magnitude. Therefore, the H, data was also increased by a 
factor of 6. The results of this normalization are shown in 
e Ct aa e 
Pagel. It 1S Seen Gaat.une H, on Argon experimental data pre- 
sented in this manner lies below the predicted curve by a factor 
of about lax 
+ 
Second, it was assumed that the experimentally measured H 
on Argon cross section was more accurate than the computed. 
Since the classical model predicts that the cross section for 
ae? 3 : + : 
H, is proportional to the square of the H cross section (Eq. oa 


+ 
the predicted curve for H, on Argon in Fig. 21 would be moved 


2 . Z : 
down by a factor of (6) or 36. gFigume 22,illustratessthig. jal 


+ 


this case, it 1S seen that the experimental measurement for H, 


on Argon lies below the cross section predicted by the classical 


model by a faectonmsof ga bouteZaie 3. 
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Figure 2l. Analysis of experimental measurement assumming a 
systematic error. 
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Figure 22. Analysis of experimental measurement assumming no 
systematic error. 
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V. CONCLUSIONS 


The experimental results given above can be analyzed in two 
distinct ways: (1) Assume that the theoretical i + Ar curve is 
correct and that a systematic error was present in the experi- 
mental measurements. (2) Assume that the experimental an + AX 
cross section 1s correct. We present below the conclusions based 
on both assumptions and our reason for believing that method (2) 
is more likely to be correct. 

If a systematic error were present and affected the measure- 
ments of both the H” and H, cross sections by the same magnitude, 
then it was seen that the classical model predicts a cross section 
for H. too large by a factor of about 14. That the actual cross 


2 


. - ‘ é 
section for H, is so much smaller than predicted 1s probably due 


to the two effects previously discussed. That is: (1) the 
Orientation probability was computed by considering the maximum 


area enclosed by. the A@ = g(Ab) curve and (2) the model neglected 


ae 
Z 


collision process. That the predicted cross section was so much 


the possibility that the H, molecule might be excited during the 


greater than observed would perhaps make the classical model 
+ : 
useful only for predicting wupper Tamats of H, cross sections. 


However in view of the agreement of our af on Argon data with 
Magnuson et.al. non. this is considered the less likely of the 
two conclusions. 

If both experimental measurements are correct, then it is seen 


that the classical model predicts an H5 cross section too large 


bg 


by a factor of about 2 to 3. This would make the classical moeaqel 
useful for predicting H, cross sections on an absolute basis. 
This is considered to be the more likely of the two conclusions. 
Assuming that the classical model is essentially correct the 
following simple refinement is proposed: Recalling equation (3) 


it was seen that the orientation probability could be written as 





where S was the area enclosed by the A@ = g(Ab) curve. Since § 
could not be determined exactly, an upper limit was obtained by 
considering the maximum values for A@ and Ab respectively when 
the other variable was zero. This in fact gave four points on 


the AB’ = g(Ab) curve, n&mely: 


Ag 


=O, ab = + Ab 
Ney = oe = - Ab. x 
AR rome = it rar 
Ab = 0, AY = - ND aes, 


A first Approximatrom tothe ay-= g(Ab) curve cane thererore we 
obtained by simply connecting these points with straight lines. 


The following sketch illustrates this. 
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The area enclosed by this first approximation is 


s=3{2(db ag = oll2p(ag_ ab = O] = & S 


max 
which leads to P = & P 
max 
: 
and 0 (0) = a ioe) a (9) 
H5 Th Ey H 


Equation (9) predicts the experimentally observed cross section 


ap : 
EOE H,, to within the accuracy of the experiment. 
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1. Computer Program for 0 


APPENDIX A 


H H5 


,(®) anda |(0). 


This computer program solves the following set of equations 


for classical differential scattemmanpin the Center of mass 


coordinate system. 


where @ 


with [ 


V(x) 





ae oD 
SinO0 dQ 
pes ol! 
Lee} 
1 |, _.v¥@) . 
2 apie 
r : oo 
O 


largest positive root 


a! 


oO 


reduced mass of the system 


velocity of incident particle in the laboratory 


system 


scattering potential 


In this case the scattering potential as given by Byatt (3) 


for Argon converted to electron voits is V(R) 


with R in units of Bohr radii (0.52928) where 


A= 411.4 
Be==-11 7 > 
C = 195.9 
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= =( gee pe 2C8 ) 
= s<21eo76 
bea = Sole / 
You- 9.625 


An outline of this program is as follows: 
1. Read in the desired incident energies, mass of proton, mass 
of target atom, disassociation energy, and potential parameters. 
2. Select an impact parameter, b. From an initial value, b is 
increased in units of .05 ave where ay is the first Bohr radius. 


3. Solve for the largest positive root of 


Oe b 
= a es O using the Newton-Raphson method. 
MV Ve ry 
we v(x b° “2 y 
4. Baalmate the antegra liga = ( (1 -YGl. > ] a) ar 
“gL x r 
if O 


O 


using the Gauss-Mahler quadrature formula 


tah Sit IT TT 
AS (cos =) 9 ( <a eE(Gocm =. )1G\(ecap meee Sot) 
al 
V(r /x; ) 2 % 
Vie TeX. ee ee 
1 Ly 2 (x _/x " 
=~ © on 1 
with x. = cos le es COS i. > a — au This procedure was taken 
. i: o> ” 22°" 7p Besa 


from Bates, Cook, and Smith ey 


L do) _ du du. 
5. From Q = 1 - 2bU, evaluate |S- |= Jap + 2b- |. a is 

obtained from the analytic expressions 
V(r) be 

U = f(x 0b) andl > => = aoe O in a straight-forward 
p/2 Vv to 

dr 
du _ 0U , QU Oo 





manner using the chain-rule of differentiating 3p —— Ele 5 se 


7x0 aon 


Ox Solve for oe = EG 
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8D 





7. Solve for the orientatpomsemebability see Oo (9) 


8. Solve for 0 


a ea H 


9. Write out incident energy, b, r_, 0, 0 POO) a P , WIE). 


Oo 


Computations are done in-units of Bohr radii, but 0 mee) and 


H» 


a 
o (9) are converted to cm at the end of the program. The 


Hy 


following is a list of key symbols used in the computer program: 


DENER 
RHO(J) 


EV(K) 


~ 


ET 


DV 


DEG(J) 


SIGMA 


DisassociatimEnergy, D, in units of eV. 
Impact parameter, b, in units of Bohr radii. 
Incident ion energy in units of eV. 
5 BYU) Getinition 
yt 

V(R) R : : : 

AT convenient form of the potential in 
this program 


av 
dR 


scattering angle in degrees converted to 


laboratory system 
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2. Calculation of Ion Trajectory and Solvay...) 

As discussed in Section III E, the scattered ion trajectories 
and solid angles must be computed numerically. Since much of the 
computer program that accomplished this is peculiar to the Naval 
Postgraduate School IBM 360-67 computer installation, only an 
outline of the program is presented here. 

Recall that the axial component of the magnetic field on the 
magnetic field axis was measured and expressed as a polynomial 


of the form i 3 


Bush [4] shows that a vector potential for this axially symmetric, 


non-uniform magnetic field can be written using this expression 


for Bi(z, =O ) as: 
oo n 2nt+1 we zr = 0) 
AS, 2 =a es = aad 
g\ 4 i n=@ n! (nie (a en 
Oz 
A (© 5z) = 50 
A, (452) = O 


Using this vector potential in the Lorentz force equation 


al 


F=4q x B=qv X(V x A)S 


it follows thats 











2 A A, 

_ eae ag ag °F Ag 

oe yee ($2 q 3 Ox | a) 
2 OA el EA 

: dr\ /ad\ | ag | _ dz\/°"g)\ 1 a PD} (2) 

"Soi 2 dita as “$I ane (EV a Je Ge Woke 
2 OA 

= d ag 

Fi = m 2 = q ao Aa (3) 
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| a 4g | | 
Equation (2) can be reduced to ee Using this, 


equations (1) and (3) can be combined to eliminate the time 


dependence to 





ace 

ie 2 3 

= oe 2 _ Hip hes _ aA 
fae (S2y*} | ar? Ea 


This is a second-order, non-linear differential equation whose 
solution describes the trajectory of the scattered ion in the 
magnetic field: In rder™to integweate it Wimeriecalm,, it 1s 
expressed as two first-order, simultaneous differential equations 
as follows. 


Let Vy = ir 





uU 
N 


O 


a U 
> 
aS 


Then equation (4) may be equivalently written as 


ele 

dz "2 (5) 
dy 

Za 3 2 iz 6 
Seo a a ) 
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Equations (5) and (6) were integrated by the DHPGC subroutine 
which  . a Hamming's modified ee ee method. 
The solution A (5) and (6) yields the trajectory and hence Zo 
and Oo: Now in a Single computer run, the trajectory equations 
were solved a number of times for various scattering angles. 
Typically this would be for @ from 35° to 55° in 1° increments. 
A second subroutine was then called on to generate a seventh 
order least squares polynomial fit of 9 = £(Z,)- From On and 
the. detector geometry, AZ was computed. The solid angle could 


then be computed by the expression 


~ 290\ _. 30 
GANZ =a 303) sin 0 lsz_ | AZ. 


The computer output then tabulated 9, Z , AZ, 9 , and aki as 


well as a plot of the trajectory (see Fig. 12) if desired. 


7O 
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